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INTRODUCTION 

A symposium of this type generally attempts to collect two types of papers: papers that 
present new work and those that try to summarize and act as pathfinders. The latter class of 
paper is generally written by harried researchers who no longer have time to verify that their 
ideas seldom work as initially conceived. The reader can readily surmise that this paper belongs 
to the latter category. 

The first part of this paper will simply discuss some considerations toward developing nu- 
merical procedures for simulating viscous compressible flows. Both Navier-Stokes and boundary- 
layer field methods are considered. For the most part, this discussion draws on our own ex- 
periences using finite difference procedures and does not pretend to cover the literature. The 
remaining parts of this paper will deal with two topics with which the authors are involved: 
1) a simple formulation for the three-dimensional boundary-layer equations in arbitrary gener- 
alized coordinates, and 2) discussion of a technique which we have referred to as the fortified 
Navier-Stokes (FNS) approach. Both topics are directed towards a long-term effort to build 
a general purpose compressible flow solver that can optionally take advantage of approximate 
solution methods both to improve accuracy and efficiency. The three-dimensional boundary- 
layer procedure in generalized coordinates fulfills a need that we feel has not been adequately 
satisfied, that is, a simple boundary-layer procedure that can share grids, turbulence models, 
and even variables with a Navier-Stokes procedure (specifically, a boundary-layer code that 
does not require special assumptions such as coordinate orthogonality in any one direction). 
The fortified Navier-Stokes scheme is a procedure that allows approximate techniques or any 
known information to be incorporated into the Navier-Stokes scheme in a rather benign way. 
If the approximations break down, the method can still proceed as a less efficient conventional 
Navier-Stokes solver. 

Throughout this paper the question of turbulence modeling will be avoided. Although it 
is generally acknowledged that current turbulence models are inadequate for complex flow field 
simulation, our own intuition is that other problems must be dealt with first. Too often, poor 
results have been blamed on turbulence models when error in specifying boundary conditions 
and poor grid resolution are the real culprits. Certainly the problem of turbulence modeling 
cannot be properly assessed until these errors are eliminated and we can routinely deal with 
complex three-dimensional flows. 
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BACKGROUND 


Numerical algorithms for viscous flow simulation are designed with several goals in mind. 
The simulation must be accurate, efficient, and simple and robust so as to minimize the amount 
of engineering manhours time needed to obtain a solution. Because we do not have sufficient 
computer resources (and perhaps sufficiently well conditioned algorithms) compromises have to 
be made. Viscous flow simulation is often carried out using either a coupled boundary-layer 
inviscid-flow model or the Reynolds- averaged Navier-Stokes equations. For supersonic flow 
and channel flows, parabolized Navier-Stokes equations which can be marched in attached flow 
regions have been used as well. 

Viscous-inviscid interaction schemes for airfoil analysis have been one and even two orders of 
magnitude faster than Navier-Stokes schemes, but they have proved to be difficult to generalize 
to complex three-dimensional separated flow situations and require considerable amount of 
engineering time to set up. Moreover, the limitations of the boundary-layer equations are not 
fully understood. For this reason, Navier-Stokes schemes are generally used to simulate complex 
flows even though they require considerably more computer work per point and usually lead 
to stiffer sets of equations to solve. Consequently, in viscous flow simulation using the Navier- 
Stokes equations the computational fluid dynamicist may spend the majority of his or her time 
trying to make the numerical algorithm more efficient. In viscous flow simulation using matched 
boundary-layer and inviscid flow equations the main task is to make the codes general enough to 
reproduce complex flow yet simple enough to economize on the engineering resource. Thus one 
approach entails making a general scheme more efficient, the other entails making an efficient 
scheme more general. 


Navier-Stokes 

The Navier-Stokes equations are considered to be adequate for most aerodynamic applica- 
tions and can be used for regions with massive flow separation. Navier-Stokes algorithms that 
are capable of treating complex geometries use finite difference methods or related finite volume 
and finite element methods. Because a limited number of grid points are available, regions of 
massive separation are generally computed in high Reynolds number flow as rotational inviscid 
flow even though the viscous terms may be retained throughout. Compared to other meth- 
ods of simulating strongly interacting flows, Navier-Stokes simulation is straightforward and 
does not require as extensive engineering ‘set-up’ time as viscous-inviscid interaction methods. 
However, Navier-Stokes schemes tend to be computer time and storage intensive. As a result, 
considerable effort is being expended to improve the computational efficiency of Navier-Stokes 
codes. This effort has, in fact, been successful to the point that algorithm improvements have 
been keeping pace with improvements in computer hardware. Overall computational efficiency 
has been improved in a number of ways including use of space varying time steps relaxation 
variables, use of a sequence of mesh refinements to establish approximate solutions, inclusion 
of various ways to reduce inversion work, improved numerical dissipation terms, more accurate 
difference approximations, and more implicit treatment of terms. References 1 to 21 provide 
one or two pointers to some of the commonly used Navier-Stokes code and descriptions of some 
the algorithmic changes alluded to above. 

Two examples of complex three-dimensional flow simulations carried out with a Navier- 
Stokes scheme are shown in figures 1 to 4. Figures 1 to 3 show results for jet in a cross flow 
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which were extracted from reference 22, while figure 4 shows preliminary results (‘oil-flow’) for 
supersonic flow over the integrated Space Shuttle vehicle (unpublished work of Buning, Chiu, 
Obayashi, Rizk, and Steger of a NASA Ames Research Center Space Shuttle flow simulation 
group). Both results were carried out on one processor of a CRAY 2 computer, used about 
a quarter of a million grid points, and took from 4 to 8 hours of computer time. Although 
the jet in cross flow represents a simple configuration, the flow field is complex and more grid 
resolution (or better numerics) is needed to better resolve the flow features which are unsteady. 
The simulation of the integrated shuttle vehicle uses a chimera composite grid approach (ref. 
23), but binding linkages, plume effects, external fuel feed lines etc. have not yet been modeled. 

Viscous-Inviscid Interaction 

Solution algorithms using viscous-inviscid interaction schemes are perhaps equally divided 
between those which use integral boundary-layer methods and those that use either finite dif- 
ference or finite element methods — field boundary-layer methods. Integral boundary-layer 
methods are quite fast and with use of good correlation functions, they can give good accuracy. 
However, they require modeling of the flow physics, in addition to the turbulence model. As 
a result it is difficult to decide whether errors are due to turbulence modeling or the integral 
method modeling. The integral methods also require extensive correlation with experimental 
data. For these reasons we have preferred field methods over integral methods. (It has also been 
our experience that once the the boundary-layer and inviscid flow equations are coupled, the 
computer time savings accrued by using an integral boundary-layer scheme over a field scheme 
is minimal.) 

Compared to Navier-Stokes schemes, viscous-inviscid interaction methods require signif- 
icantly less computer time to simulate viscous flow with up to small regions of separation. 
For high Reynolds number flow, accuracy can be equal to Navier-Stokes simulation, however, 
depending on the problem, engineering ‘set-up’ time can be higher. For strong viscous inter- 
action, current boundary-layer based schemes tend to break down, although we believe that if 
the boundary-layer equations can be properly interacted with the inviscid rotational flow equa- 
tions, for many problems they can give high Reynolds number flow results comparable to the 
Navier-Stokes equations. 

Viscous flow simulation using boundary-layer procedures requires special adaptations and 
thus has an intensive user workload. For example, if separation is encountered, the boundary- 
layer algorithm must be converted from the direct to the inverse mode to avoid singular behavior 
at the separation point. A scheme to update the inverse mode quantity, generally r w or 6 *, must 
then be provided so that the pressure predicted by the boundary-layer and inviscid schemes are 
consistent. In conventional viscous - inviscid interaction, the influence of the viscous layer 
is imposed by displacement thickness or transpiration velocity. An advantage of using either 
transpiration or an effective body displacement is that the inviscid flow can be computed using 
a grid that is relatively coarse compared to what is needed for the viscous grid. Moreover, the 
inviscid grid need not necessarily be body conforming in the way that a high Reynolds number 
Navier-Stokes grid must be. However, matching to a highly rotational inviscid outer flow is not 
straightforward; and, in many cases, it is necessary to account for the viscous flow curvature 
and the pressure gradients of the shear layers. 

As an example of the inviscid-viscous interaction approach, the results of computing flow 
about the RAE 2822 airfoil at M <*> = 0.73 and Ci = 0.803 are shown. Figures 5 - 7 show 
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comparisons of the C P ,C /, and 8* distributions found experimentally (ref. 24), those computed 
by the interaction method (ref. 25), and those computed by Mehta (ref. 26) using a Navier- 
Stokes code. Mehta performed his computations at a = 2.79° (and computed a Ci = 0.793); 
the present computations were performed at a = 2.81° to match the measured Ci — 0.803. The 
present results are in good agreement with both the Navier-Stokes and experimental results. 
These simulations were obtained on 223x31 inviscid and 223x50 viscous grids using a quite fast 
transonic full potential code (ref. 27) and equally fast boundary-layer finite difference code 
(ref. 28). The calculations were carried out in 1984 and for the cases presented, the required 
Cray-XMP CPU time was 7 to 15 sec, and on the average, 0.0006 sec/grid point were required 
to obtain a converged solution. In contrast, our most efficient thin-layer Navier-Stokes code at 
that time (ref. 7) required about 0.03 to 0.06 sec/grid point to obtain a converged solution. 

Zonal Equation Methods 

The traditional choice in viscous flow simulation is between either a viscous-inviscid in- 
teraction approach using simplified equations or a Navier-Stokes scheme used throughout the 
field. In practise, however, various kinds of zonal methods are often used. A zonal equation 
method generally entails regional use of simplified methods with a Navier-Stokes scheme so as 
to speed up iterative convergence, reduce computational work, or even improve the numerical 
accuracy of the Navier-Stokes simulations. Simple zonal ideas are almost always used in high 
Reynolds number flow Navier-Stokes simulations. For example, for high Reynolds number flow, 
the Navier-Stokes equations revert to Euler equations in any region in which sufficient grid res- 
olution is not provided, so an obvious zonal scheme is to neglect viscous terms in these regions. 
Special treatment at wall boundaries using ‘wall functions’ (ref. 29) have also been employed. 

In a zonal approach the simplified equations must be able to readily share variables and 
mesh points with the global scheme. Otherwise, sometimes awkward conversions of variables or 
messy interpolations between equation sets are required. 

One kind of zonal approach is to simply embed Navier-Stokes equations in complex (usually 
separated or highly vortical) regions of the flow and use simplified equations elsewhere. This 
requires some care in establishing zones, and the zonal interfacing can be tedious. Alternately 
a zonal approach can use a global Navier-Stokes and Euler formulation throughout, but use 
simplified equations embedded into the global scheme to reduce computational work or improve 
accuracy. The subtle difference between this zonal approach and the former (which is perhaps 
chiefly observed in the coding) is that, if the approximation behind the simplified equation is no 
longer valid, simulation can still proceed with the more general scheme. The fortified Navier- 
Stokes scheme (to be highlighted in a latter Section) and related embedding schemes which use 
vorticity, potential, and other like forms of the equations (refs. 30-32) are pronounced examples 
of the latter kind of zonal scheme. An advantage of a general zonal approach is that it is possible 
to directly compare the effects of using the simplified equation set to the general equation set 
and thus learn the limitations of the approximation. As experience is gained, it is then likely 
that the simplified equation sets will be able to assume more of the computational work. 
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GENERALIZED BOUNDARY-LAYER EQUATIONS 

The boundary-layer equations are restricted to high Reynolds number flow and require 
the use of a body conforming coordinate system. In developing the usual boundary-layer 
equations, both the independent variables and the dependent velocity variables are transformed 
to the new body conforming coordinates. For body surfaces with little curvature, the boundary- 
layer equations cast in terms of the new dependent variables essentially simplify back to a 
flat plate or Cartesian-like form of the equations along a developed surface. If the body has 
appreciable curvature, however, the equations require additional terms and can be difficult to 
solve numerically. They are particularly more complex if a nonorthogonal coordinate system is 
used, yet for many applications it is difficult to generate an orthogonal coordinate system along 
the body surface. 

The thin-layer Navier- Stokes equations use body conforming coordinates and only viscous 
terms in the normal-like direction are retained. With the pressure specified in the viscous layer, 
the momentum equations of the thin-layer Navier- Stokes equations are readily solved for the 
original Cartesian velocity components provided only the independent variables are transformed 
to body coordinates. This is because with pressure specified, the thin-layer momentum equa- 
tions are essentially equivalent to an uncoupled set of scalar convection-diffusion equations. 
Moreover, the boundary-layer approximation of specifying the pressure in the viscous layer is 
equivalent to solving a linear combination of the thin-layer momentum equations - the combi- 
nation needed to transform the Cartesian velocity components u,v, and w into an equivalent 
normal-like velocity (see ref. 33 ). Consequently, although three thin-layer momentum equa- 
tions in transformed independent variables are readily solved for u,v, and w , for pressure pre- 
scribed in the viscous layer only two linear combinations of the u,u, and w velocity components 
should be retained - the combination that is equivalent to the velocities along the body surface 
coordinates. Thus equivalent boundary-layer equations can be formulated from the thin-layer 
equations for pressure specified in the viscous layer; and, because only the independent variables 
are transformed, complex coordinate source terms can be avoided. Moreover, this alternate form 
of the boundary-layer equations does not require that any of the coordinates be orthogonal, and 
software (grids, boundary condition routines, etc.) developed for many Navier-Stokes schemes 
can be readily utilized. 


Governing Equations 

A form of the boundary-layer equations for general curvilinear coordinates £ = £(x,y, z), 
77 = rj(x,y,z ), and £ = £(a;,y,z) is given by reference 33 (with 77 chosen away from the surface): 


77 Momentum 

V77 • (V£p € + Vt7p^ + V<p<) =pU(ud i fj x + vd^fj y + wd(fj z ) 

+ pV(ud v fj x + vd v f) y + wd v fj z ) 
+ pW{ud(fj x + vd(fj y + wd(fj z ) 

or for very dominant V77 • V77 


Pv — 0 


(la) 


( 16 ) 
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£ and ( Momentum, 

Predict values of q = (u,v,wY from 


pu t + pUu i + pVu v + pWu i +{£ x p ( +VxP v + CxPc) = J5 n [J _1 (/im iu v + (p/Z)m 2 rj x )] (2a) 

pv t + pUvt +pVv v + pWv< + (t y pt +VyPv +CyP<) = JdnlJ^ifJtmiVr, + ((i/3)m 2 T}y)] (26) 

pw t + pUw{ + pVw v + p\Vw( + (Z z pt + r) z p v + ( z p<;) = J 4- (p/Z)m 2 ri z )} (2c) 
and retain the two linear combinations 

U = V(-q 

W = V C-q 

Energy 


pH t + puH{ + pvH v + pwH( - p t = Jd v [J 1 (pm 1 m 3 + (p/3)m 2 (Tj x u + rj y v + -q z w))} (3) 


where 


Constitutive 


where 


Continuity 


m 1 =r)l+T] 2 y + r] 2 z 

m 2 = TJzUr, + TJyVr, 4" J|zW, 

m 3 = [( u 2 4- v 2 4- w 2 ) r) ]/2 4- Pr -1 ( 7 - l) -1 (a 2 )^ 


P _ PT 00 
Poo P 00 T 


T_ _ ( 7 - l) f g _ (u 2 4-v 2 4-W 2 ) 

P -2 !■““ O J 


(4a) 


(46) 


( J“V)r + + (J-'PV)* + (J-'pWh = 0 (5) 

Throughout p, p, and T denote fluid density, pressure, and temperature, u,v,w are Cartesian 
velocity components referenced to an inertia system, H is the stagnation enthalpy, a is the sound 
speed, and p is the viscosity coefficient. The transform Jacobian is given by J = 1 |, and 


552 



the overbar on rj x , rj y , and r] z denotes scaling, r) x — r\ x j + rf% etc. The operator V 

denotes the vector gradient so that V77 • q = rj x u + rj y v -f r] z w — V, and V is the vertical 
contravariant velocity. 

The six equations (77 momentum (1); two linear combinations of (2a), (2b), and (2c) to 
form U and W\ energy (3); state (4a); and continuity (5)) can be used to determine the six 
variables p, U, W, H, p, V. The Cartesian velocity components are then obtained from 


u 

v 

w 


x i 

X V 

X C 

Vi 

Vv 

y< 

.H 

z v 

z c. 



( 6 ) 


Numerical Scheme 

This formulation has been applied to steady state problems using the time-like boundary- 
layer scheme reported by Van Dalsem and Steger in reference 34. For a prescribed edge pressure, 
the equations are solved iteratively from an assumed initial state in the following way with 
p v = 0. Using central spatial differencing in 77 and upwind differencing in £ and £, equations 
(2) to (3) are used to update u,v,w , and H. As pressure was already obtained using one 
linear combination of the momentum equations, only two linear combinations of the momentum 
equations can be used to determine the velocities. Thus, U and W are formed from u,v, and w. 
The equation of state, (4a), is used to update p , with T defined from equation (4b). The third 
linear combination of u,v,w is obtained by solving continuity for V using already updated U, 
W, and p. The continuity equation has central space differencing in £ and £, and trapezoidal 
differencing in 77 so that V is obtained by marching outwards from the body surface. The final 
updated form of the Cartesian velocity components is then obtained from updated U, V, W 
using equation (6). This process is then repeated until convergence. 

The boundary-layer equations are solved with either specified initial data profiles or the 
stagnation region can be captured with the upwind differencing as part of the iteration process. 
At the outer edge in 77, variables can be specified if known, or, for specified edge pressure, the 
edge variables can be computed as part of the boundary-layer method. In this case the pVd v 
terms in momentum and energy are backward differenced away from the edge if V > 0, or 
dropped if V <0. Likewise the viscous terms must be dropped at just the outer edge point. 
Periodic, symmetry, or outflow conditions are used at all other £ and £ boundaries at which 
initial data are not specified. 


Numerical Results 

Boundary-layer calculations to verify this formulation have been carried out on a 6 : 1 
prolate spheroid at 10° and 30° angle of attack and at two different Reynolds numbers, 7.7 x 10 6 
and 43 x 10 6 based on diameter. The geometry and a top view of the computed limiting 
streamlines for the 30° angle of attack case are shown in figure 8. Streamwise separation was 
avoided by confining the boundary-layer calculation from x/L = 0.02 to x/L = 0.8 where x is 
the axial distance and L is the body length. The initial profiles in £, edge conditions, turbulent 
eddy viscosity values, and even the grid itself were taken from the Navier-Stokes calculations 
described in reference 35. Throughout, the simplified 77 momentum equation, p v = 0, was used. 


553 



As reported in references 33 and 35, good agreement between the boundary-layer method, 
the Navier-Stokes calculation and even experiment is obtained for attached flow regions. For 
latter reference, a boundary-layer computed profile at x/L = 0.65 and <fi = 120° ( <j> = 0° 
is the windward plane of symmetry) is shown in figure 9 along with the computed thin-layer 
Navier-Stokes result. 

Calculations were also carried out for the 10 ° case in which the boundary-layer code was run 
from x/L — 0.8 to x/L = 1.1 - the latter x/L location is back on the supporting sting which was 
modeled in the Navier-Stokes simulation. In this region streamwise separation occurs. As shown 
by the separated profile at x/L = 0.96and <f> = 120°, figure 10, the boundary-layer equations 
solved in the direct mode break down very rapidly unless the outer ‘edge’ is taken very close 
to the wall. To prevent breakdown the constant height edge had to be taken so sufficiently 
close that for this profile the flow is not just rotational but still fully viscous. Breakdown of the 
boundary-layer equations in direct mode (i.e. pressure specified) is not surprising. Of interest, 
though, is that the boundary-layer equations do give good results in the direct mode provided 
that good outer edge values (not irrotational values) are given sufficiently close to the wall. 

FORTIFIED NAVIER-STOKES 

In a typical high- Reynolds number Navier-Stokes simulation the fine-grid resolution is gen- 
erally provided in a thin zone near the body surface, and the outer flow is effectively resolved 
as rotational inviscid flow (unless the turbulent coefficients are quite large, in which case exten- 
sive modeling is required). Because of the fine-grid resolution required near the body, a given 
algorithm often operates much less efficiently on the Navier-Stokes equations than it does on 
the Euler equations, even though the work per step may be similar and the viscous terms may 
enhance stability. However, on this same refined viscous grid, the boundary-layer equations can 
be efficiently and accurately solved. As a result one can speculate that by using a boundary- 
layer algorithm near the wall, it may be possible to significantly improve the productivity of the 
Navier-Stokes algorithm. A convenient way of imposing the boundary-layer equations in this 
way is by a general method which we have termed fortified Navier-Stokes (FNS). 

Fortified Approach 

As an alternative to the traditional “zonal” concept of solving the various flow zones on 
separate grids and patching the zones together, the Navier-Stokes equations can be applied 
throughout and the solution of simplified or subset equations can be embedded into the more 
general or global algorithm by means of a forcing function (fig. 11). In this way, the approximate 
solution scheme is used to enhance or “fortify” the more general scheme so as to improve 
the overall accuracy or efficiency. For example, the thin-layer Navier-Stokes equations can be 
represented in a form 


dr{Q) + di(F) + d v (G ) + d^H) = Re~ A d-qS 

where Q is the solution vector, F, G , and H are flux terms, and S is the viscous stress term in 
77 . These equations can be fortified through the addition of a dissipative forcing term as (refs. 

36-38) 
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dr(Q) + d ( (F) + d v (G ) + d^H) = Re~ 1 d v S + C(Q bl - Q) (7) 

where Q&i denotes Q formed from the boundary-layer solution. The positive definite operator 
C is chosen such that — CQ is a damping term and is easily invertable, for example 

c = xl (8a) 

with I the identity, or 

C — + driv + (86) 

The parameter x is picked as a large positive coefficient that blends to zero outside of the forced 
region of the viscous layer. Inside the viscous layer Q b i approaches Q so no error is generated 
by the addition of the term C(Qbi — Q )• For a frozen QbU such an implicitly imposed forcing 
term can be used to favorably alter the stability properties of numerical solution algorithms. 
(Analysis of this approach has been presented in references 37 and 38.) Outside of the forced 
region of the viscous layer x is zero and only the thin-layer Navier-Stokes equations are solved. 

By adding subset equation forcing terms to the Navier-Stokes equations, the Navier-Stokes 
and subset equations interact strongly over entire regions rather than just at interface bound- 
aries. Moreover, the subset equations can be applied selectively to only those regions where 
they are clearly valid. In principle, if during an iterative solution process the approximation 
method become suspect in a certain region, the forcing can be turned off, and the region can be 
resolved with only the global Navier-Stokes scheme. Because of this flexibility, the generality of 
the Navier-Stokes equations is retained with the FNS approach, while some of the efficiency of 
the subset algorithms is recovered. 

Numerical Implementation 

The fortified terms have been implemented in an implicit approximately factored algorithm 
for the three-dimensional thin-layer Navier-Stokes equations written in general coordinates. The 
basic two-factor solution algorithm (refs. 14, 15) uses central spatial differencing in the T) and £ 
directions and upwinding in the £ direction. The forcing term with C = x^ has been added to 
the basic algorithm as (refs. 37,38) 

7(1 + h x ) + M£(i+) n + h6 c t 7 n - hRe- 1 6 < J~ 1 M n J - £>;| c ] 
x [1(1 + hxT 1 x [l(l + h x ) + hS((A~) n + h6 v B n - DiU] AQ n = (9) 

- At{^[(F+r - F+ ] + */[( F") n - F" ] + *„(G n - Goo) + *<( H" - H^) 

- Re-'Sd S n - So,)} - (D e |, + D e | C )(Q» - Qoo) + h X ( Q/ - Q n ) 

Here h = At (first order in time), or /i = y, (second order in time) and a free stream base 
solution is subtracted out to improved accuracy in the far field. The operators and 6 are 

backward and forward three-point difference operators. The flux F has been eigensplit and the 
matrices A,B,C, and M result from local linearization of the fluxes about the previous time 
level. Because central-space-difference operators are used in 7] and £, implicit Di and explicit 
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D e numerical dissipation terms are included in equation (9). With this implementation, x 
can become very large without concern for large factorization errors, and x contributes to the 
diagonal dominance of both left-hand factors. 

Results 

The FNS method has been tested on a simple geometry which roughly models an infinitely 
swept wing - a 7.75% thick sine-wave bump-on-a-wall with a 35° leading edge sweep. The flow 
was computed with a thin-layer Navier-Stokes algorithm (refs. 14,15) and with the FNS scheme 
using the Van Dalsem and Steger boundary-layer field method (ref. 34). The flow was first 
computed using the standard Navier-Stokes algorithm alone on both a fine mesh (29 points in 
the flow direction, 50 points normal to the wall, 5 points in the span direction) and a coarse mesh 
which has only 20 points in the critical normal direction. The same minimum normal spacing 
at the wall was used in both the coarse- and fine-mesh computations. This flow was also 
computed with the FNS approach using the coarse Navier-Stokes mesh (20 points in the normal 
direction) and, near the wall, a superimposed fine boundary-layer mesh (50 points in the normal 
direction). In all the FNS computations presented here, % was proportional to the vorticity; 
hence, it is automatically large within the boundary-layer and rapidly drops to zero near the 
edge of the boundary-layer. The drag history in figure 12 shows that the coarse-grid standard 
Navier-Stokes computation does not predict the drag accurately, and that the FNS method 
obtains essentially the same drag level in 50 iterations that the fine-grid standard Navier-Stokes 
computation reached in 400 iterations. The computed near-surface particle traces are shown 
in figures 13a-c. Both the FNS method and the fine-grid standard Navier-Stokes computations 
predict a constant chord-line separation line, whereas the coarse-grid standard Navier-Stokes 
computation does not quite capture this qualitative feature. 

The swept infinite wing geometry described above was modified to yield a three-dimensional 
flow by reducing the aspect ratio of the wing to one. Also, to resolve spanwise changes, the 
grid dimension in this direction was increased from 5 to 15. The resulting particle traces (figs. 
14a-d) show the same trends as described for the infinite span example. Also shown in figure 
14d is the result obtained when the viscous terms and no-slip boundary conditions are removed 
from the global numerical algorithm (thus making it an Euler formulation). In this case, the 
entire influence of viscosity must be carried by the boundary-layer algorithm, which is not a 
difficulty for this case. The drag history versus CRAY-XMP CPU time for these computations 
is presented in figure 15. As before, the coarse-grid standard Navier-Stokes computation is not 
accurate, and the fine-grid standard Navier-Stokes computation is expensive, whereas the FNS 
(and Euler) computations yield the same drag level as does the fine-grid standard Navier-Stokes 
computation, but for one-twentieth of the cost. 
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CONCLUDING REMARKS 


Because efficient viscous-inviscid interaction methods have been difficult to extend to com- 
plex three-dimensional flow simulations, Navier-Stokes procedures are more frequently being 
utilized even though they require considerably more work per grid point. It would seem a mis- 
take, however, not to make use of the more efficient approximate methods in those regions in 
which they are clearly valid. Ideally we should like a general purpose compressible flow solver 
that can optionally take advantage of approximate solution methods both to improve accuracy 
and efficiency. Some potentially useful steps toward this goal have been described in this paper 
- a generalized three-dimensional boundary-layer formulation and the fortified Navier-Stokes 
procedure. Further work remains in merging these steps and more. 
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a. (laminar) 



b. (turbulent) 


Figure 1. Jet in ground effect issuing into a crossflow with y 33 - = 0.223 and a nozzle height 
located 3 jet diameters above the wall. 
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Figure 2. Centerline C p distributions for the jet in ground effect with y 33 ^ = 0.223 and a nozzle 
height located 3 jet diameters above the wall. 
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Figure 3. Time history of particle traces indicating the interaction of the ground vortex and 
the ring vortex. 
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Figure 4. Simulated oil flow (limiting streamlines) compared to experiment for the integrated 
space shuttle at M oo — 2.0, a - —4°. 
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Figure 5. Computed and experimental pressure distributions for RAE 2822 airfoil at M ^ = 
0.73, R eoo = 6.50 x 10 6 , Ci = 0.803. 



Figure 6. Computed and experimental C/\ e distributions for RAE 2822 airfoil at M a 0 = 
0.73, = 6.50 x 10 6 , Ci = 0.803. 
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Figure 7. Computed and experimental 8* distributions for RAE 2822 airfoil at M 0 0 = 0.73, Re oo = 
6.50 x 10 6 , <7, = 0.803. 




Figure 8. Top view of 6 : 1 prolate spheroid showing thin-layer Navier-Stokes and boundary-layer 
(0.02 < x/L < 0.80) computed limiting streamlines for M 0 0 = 0.25, a = 30°, and Re = 43 X 10 6 . 
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Figure 9. A typical computed attached profile (contravariant U versus normal height) at x/L = 
0.65 and </> = 120° on the prolate spheroid at 10° angle of attack using boundary-layer and 
thin-layer Navier-Stokes equations. 
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Figure 10. A separated profile at x/L = 0.96 and <f> — 120° on the prolate spheroid at 10° angle 
of attack showing breakdown of the direct mode boundary-layer computation as the outer edge 
is moved away from the wall. 
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Figure 13. Computed near-surface particle traces for the bump flow: a) Fortified Navier-Stokes; 
b) standard Navier-Stokes (fine mesh); c) standard Navier-Stokes (coarse mesh). 
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Figure 14. Computed near-surface particle traces for the turbulent flow over a 7.75% sine wave 
bump with a 35° leading-edge sweep, Re = 5 x 10 s , = 0.5, and AR=1: a) Fortified Navier- 

Stokes; b) standard Navier-Stokes (fine mesh); c) standard Navier-Stokes (coarse mesh); d) 
Fortified Euler. 
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Figure 15. Drag coefficient history versus CRAY-XMP processor time for flow of figure 14 
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